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Abstract 

Using our present knowledge on effective hadronic theories, short-distance 
QCD information, the 1/Nc expansion, analyticity and unitarity, we derive 
an expression for the pion form factor, in terms of m^, mx, Mp and f-,^. 
This parameter-free prediction provides a surprisingly good description of 
the experimental data up to energies of the order of 1 GeV. 
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1. Introduction 



The QCD currents are a basic ingredient of the electromagnetic (vector) and 
weak (vector, axial, scalar, pseudoscalar) interactions. A good understanding of 
their associated hadronic matrix elements is then required to control the impor- 
tant interplay of QCD in electroweak processes. Given our poor knowledge of the 
QCD dynamics at low energies, one needs to resort to experimental information, 
such as e~^e~ hadrons or semileptonic decays. 

At the inclusive level, the analysis of two-point function correlators con- 
structed from the T-ordered product of two currents, has been widely used to 
get a link between the short-distance description in terms of quarks and gluons 
and the hadronic long-distance world. In this way, information on fundamental 
parameters, such as quark masses, ag or vacuum condensates, is extracted from 
the available phenomenological information on current matrix elements. The lack 
of good experimental data in a given channel translates then into unavoidable 
uncertainties on the obtained theoretical results. 

At very low energies, the chiral symmetry constraints^ are powerful enough 
to determine the hadronic matrix elements of the light quark currents.i'ii'i Un- 
fortunately, these chiral low-energy theorems only apply to the threshold region. 
At the resonance mass scale, the effective Goldstone chiral theory becomes mean- 
ingless and the use of some QCD-inspired model to obtain theoretical predictions 
seems to be unavoidable. 

Obviously, to describe the resonance region, one needs to use an effective 
theory with explicit resonance fields as degrees of freedom. Although not so 
predictive as the standard chiral Lagrangian for the pseudo- Goldstone mesons, 
the resonance chiral effective theory^ turns out to provide interesting results, once 
some short-distance dynamical QCD constraints are taken into account.! 

We would like to investigate how well the resonance region can be understood, 
using our present knowledge on effective hadronic theories, short-distance QCD 
information and other important constraints such as analyticity and unitarity. 
Present experiments are providing a rich data sample (specially on hadronic r 
decays), which can be used to test our theoretical skills in this important, but 
poorly understood, region. 

Our goal in this first letter is to study one of the simplest current matrix 
elements, the pion form factor F{s), defined (in the isospin limit) as 

{n\-\dru\iD) = V2F{s) {p^- -p^oY, (1) 

with s = = {Pt^- +p^o)^. At s > 0, F{s) is experimentally known from the 
decay u^tt^ti^ and (through an isospin rotation) from e+e" tt+tt^, while 

the elastic e~7r+ scattering provides information at s < 0. 

Theoretically, the pion form factor has been extensively investigated for many 
years. Thus, many of our results are not new. We want just to achieve a descrip- 
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tion of F{s) as simple as possible, in order to gain some understanding which 
could be used in other more complicated current matrix elements. 



2. Effective Lagrangian Results 

Near threshold, the pion form factor is well described by chiral perturbation 
theory (ChPT). At one loop, it takes the formi: 



n 967r2/2 
where the functions 
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(2) 



A{mlls, mlli,^) = In (ml/fi') + M - | + a|, In f (3) 

^ ^ s 3 \ap — 1/ 

contain the loop contributions, with the usual phase-space factor 



ap = ^l-4ml/s, (4) 

and LKfi) is an O(p^) chiral counterterm renormalized at the scale fi. The mea- 
sured pion electromagnetic radius,! (r^)^* = (0.439 ± 0.008) fm^, fixes Ll{Mp) = 
(6.9 ±0.7) X 10-3. 

A two loop calculation in the SU (2) SU (2) theory (i.e. no kaon loops) has 
been completed recently 

Using an effective chiral theory which explicitly includes the lightest octet of 
vector resonances, one can deriveQu the leading effect induced by the p resonance: 

where the couplings Fy and Gy characterize the strength of the p7 and pun 
couplings, respectively. The resonance contribution appears first at the next-to- 
leading order in the chiral expansion. This tree-level result corresponds to the 
leading term in the 1/Nc expansion, with Nc = 3 the number of QCD colours. 
Comparing it with Eq. (|^), it gives an explicit calculation of Lg in the Nc oo 
limit; chiral loops (and the associated scale dependence) being suppressed by an 
additional 1/Nc factor. 

Eq. (|^) was obtained imposing the constraint that the short-distance be- 
haviour of QCD allows at most one subtraction for the pion form factor.! How- 
ever, all empirical evidence and theoretical prejudice suggests that F{s) vanishes 
sufficiently fast for s — > oo to obey an unsubtracted dispersion relation. If this 
is the case, one gets the relationi (at leading order in and if higher-mass 

states are not considered) FyGy / = 1, which implies the well-known Vector 
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Meson Dominance (VMD) expression: 



The resuhing prediction for the 0(p'^) chiral coupling Lg, 

FvGv _ fl _ „ o 1 n-3 



2M2 2M2 



7.2 X lO"'' , (7) 



is in very good agreement with the phenomenologically extracted value. This 
shows explicitly that the p contribution is the dominant physical effect in the 
pion form factor. 

Combining Eqs. (|]) and (|^), one gets an obvious improvement of the theo- 
retical description of F{s): 

Ml 



F(s) 



p 



M2 - s 9Qtt^P 

p ff TV 



(8) 



The VMD formula provides the leading term in the 1/Nc expansion, which effec- 
tively sums an infinite number of local ChPT contributions to all orders in mo- 
menta (corresponding to the expansion of the p propagator in powers of s/M^). 
Assuming that Lg(Mp) is indeed dominated by the l/Nc result (0), the loop con- 
tributions encoded in the functions y4(mp/s, mp/M^) give the next-to-leading 
corrections in a combined expansion in powers of momenta and I/Nq- 

3. Unitarity Constraints 

The loop functions yl(mp/s, mp/M^) contain the logarithmic corrections in- 
duced by the final-state interaction of the two pseudoscalars. The strong con- 
straints imposed by analyticity and unitarity allow us to perform a resummation 
of those contributions. 

The pion form factor is an analytic function in the complex s plane, except 
for a cut along the positive real axis, starting at s = 4m2, where its imaginary 
part develops a discontinuity. For real values s < 4m2, F{s) is real. The imag- 
inary part of F{s), above threshold, corresponds to the contribution of on-shell 
intermediate states: 

ImF(s) = ImF(s)2^ + ImF(s)4. + ■ ■ ■ + ImF(s)x^ + • • • (9) 

In the elastic region (s < IGm^), Watson final-state theoremS relates the imagi- 
nary part of F{s) to the partial wave amplitude Tl for tttt scattering with angular 
momenta and isospin equal to one: 

ImF(s + ie) = a^T^F{s)* = e'^'^ sin 5lF{s)* = sin 5l \F{s) \ = tan 5} ReF(s) . 

(10) 
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Since ImF(s) is real, the phase of the pion form factor is the same as the phase Si 
of the partial wave amplitude T^. Thus, one can write a (n-subtracted) dispersion 
relation in the form: 

which has the well-known OmnesEi'i solution: 

F{s) = Qn{s) exp - / - ' . , 12 

[ vr J Ami z z — s — le] 



where 



^-Qn{s) = Ey^k^-m. (13) 

k=0 



Strictly speaking, this equation is valid only below the inelastic threshold (s < 
16m^). However, the contribution from the higher-mass intermediate states is 
suppressed by phase space. The production of a larger number of meson pairs is 
also of higher order in the chiral expansion. 
Using the lowest-order ChPT result 

(14) 

the integral in Eq. (|12D generates the one-loop function — sy4(m^/s,?/)/(967r^/^), 
up to a polynomial [in s/(4m^)] ambiguity, which depends on the number of sub- 
tractions applied^. Thus, the Omnes formula provides an exponentiation of the 
chiral logarithmic corrections. The subtraction function Qn{s) (and the polyno- 
mial ambiguity) can be partly determined by matching the Omnes result to the 
one-loop ChPT result (|^), which fixes the first two terms of its Taylor expansion; 
it remains, however, a polynomial ambiguity at higher orders. This means an in- 
determination order by order between the non-logarithmic part of the pion form 
factor which must be in Qn{s) and the one which must be in the exponential. 

The ambiguity can be resolved to a large extent, by matching the Omnes 
formula to the improved result in Eq. (^, which incorporates the effect of the p 
propagator. One gets then: 

F{s) = exp {^.M</s,ml/M'^)^ . (15) 

This expression for the pion form factor satisfies all previous low-energy con- 
straints and, moreover, has the right phase at one loop. 



" Obviously, the finite integration does not give rise to any dependence on y = ra^ / p?' 
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Eqs. (|T5]) has obvious shortcomings. We have used an O(p^) approximation 
to the phase shift 6l, which is a very poor (and even wrong) description in the 
higher side of the integration region. Nevertheless, one can always take a sufficient 
number of subtractions to emphasize numerically the low-energy region. Since 
our matching has fixed an infinite number of subtractions, the result (^) should 
give a good approximation for values of s not too large. One could go further 
and use the 0(j9'^) calculation of 6l to correct this result. While this could (and 
should) be done, it would only improve the low-energy behaviour, where Eq. (p!5|) 
provides already a rather good description. However, we are more interested in 
getting an extrapolation to be used at the resonance peak. 



4. The Rho Width 



The off-shell width of the p meson can be easily calculated, using the resonance 
chiral effective theory .i'i One gets: 



rp(s) 



MpS 
967r/2 



|^(s - 4ml) (^l + ^^(^ - 4w 



cr 



K 



Im 



Aiml/s,ml/M'') + -A(m^/s,m^/M 



(16) 



At s 



M^, TplMp) = 144 MeV, which provides a quite good approximation to 



the experimental meson width, F^^p = (150.7 ± 1.2) MeV. 

Eq. (p!6D shows that, below the KK threshold, the p width is proportional to 
the 0{p^) TTTT phase shift: T p{s) = Mp6l. Making a Dyson summation of the p self- 
energy corrections amounts to introduce the p width into the denominator of the p 
propagator, shifting the pole singularity to ^/s = Mp — p. However, expanding 
the resulting propagator in powers of s/M"^, it becomes obvious that the one-loop 
imaginary correction generated by the width contribution is exactly the same as 
the one already contained in the Omnes exponential. This suggests to shift the 
imaginary part of the loop function A{rn?p / s , m^, / M'^) from the exponential to 
the propagator. Including the small contribution from the intermediate state 
KK, one has then: 



F{s) 



M2 

p 



ReA{ml/s,ml/Ml 



+ -ReA{m\/s, m\/M^ 



(17) 



This change does not modify the result at order j9^, which still coincides with 
ChPT, but makes the phase shift pass trough 7r/2 at the mass of the p resonance. 

Expanding Eq. (|l^) in powers of momenta, one can check0 that it does a quite 
good job in generating the leading O(p^) contributions in the chiral expansion. 
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Figure 1: |F(s)p (in logarithmic scale) versus s/y/s. The continuous curve shows 
the theoretical prediction in Eq. (]T7|). 

The known coefficients for linear and quadratic logarithms of the ChPT result ,0 
are well reproduced in the chiral limit. While this is also true for Eq. (0), the 
expansion of Eq. ([T7|) gives a better description0 of the polynomial part of the 
0(/) ChPT result. 

5. Numerical Results 

We can see graphically in figure 1 how Eq. (|1^) provides a very good descrip- 
tion of the dataE3 up to rather high energies. The only parameters appearing in 
the pion form factor formula are m^r, itik, Mp and /^r, which have been set to 
their physical values. Thus, Eq. (|17D is in fact a parameter-free prediction. The 
extremely good agreement with the data is rather surprising. 

At low energies the form factor is completely dominated by the polynomial 
contribution generated by the p propagator. Nevertheless, the summation of 
chiral logarithms turn out to be crucial to get the correct normalization at the p 
peak. The exponential factor in Eq. ([T7| ) produces a 17% enhancement of 
at s = M^. 

The data at s > have been taken from e+e^ tt+tt"; so in the vicinity 
of the p peak there is a small isoescalar contamination due to the u resonance. 
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This contribution is well-known and generates a slight distortion of the p peak, 
which can easily be included in the theoretical formula.li^ The effect cannot be 
appreciated at the scale of the figure. 

We also obtain a prediction for the phase shift 5\. From Eq. ( PT[ ) we get: 

i;W=arctau{^^}. (18) 

For s « Mp, this expression reduces to the O(p^) ChPT result in Eq. (0). As 
shown in figure 2, the improvement obtained through Eq. (0) provides a quite 
good description of the experimental datalllS over a rather wide energy range. 
At large energies, the phase shift approaches the asymptotic limit 6l{s —>■ oo) = 
arctan {— ^M^/(967r/^)}. If only the elastic 27i intermediate state is included, 
^ = 1 and 6l{s oo) = 167°; taking into account the KK contribution, one gets 
C, = 3/2, which shghtly lowers the asymptotic phase shift to 6l{s oo) = 161°. 
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5. Summary 



Using our present knowledge on effective hadronic theories, short-distance 
QCD information, the 1/Nc expansion, analyticity and unitarity, we have derived 
a simple expression for the pion form factor, in terms of m^r, ttik-, Mp and f.^. 
The resulting parameter-free prediction gives a surprisingly good description of 
the experimental data up to energies of the order of 1 GeV. 

Our main result, given in Eq. (P!7[), contains two basic components. The p 
propagator provides the leading contribution in the limit of a large number of 
colours; it sums an infinite number of local terms in the low-energy chiral expan- 
sion. Chiral loop corrections, corresponding to the final-state interaction among 
the two pions, appear at the next order in the 1/Nc expansion; the Omnes expo- 
nential allows to perform a summation of these unitarity corrections, extending 
the validity domain of the original ChPT calculation. 

Requiring consistency with the Dyson summation of the p self-energy, forces 
us to shift the imaginary phase from the exponential to the p propagator. While 
this change does not modify the rigorous ChPT results at low energies, it does 
regulate the p pole and makes the resulting phase shift pass through 7r/2 at the 
mass of the p resonance. 

As shown in figures 1 and 2, the experimental pion form factor is well repro- 
duced, both in modulus and phase. Although the p contribution is the dominant 
physical effect, the Omnes summation of chiral logarithms turns out to be crucial 
to get the correct normalization at the p peak. 

Many detailed studies of the pion form factor have been already performed 
previouslyJ§§lli§'§iffl The different ingredients we have used can in fact be 
found in the existing literature on the subject. However, it is only when one 
combines together all those physical informations that such a simple description 
of F{s) emerges. 

Our approach can be extended in different ways (including two-loop ChPT 
results, p' contribution, . . . ). Moreover, one should investigate whether it can be 
applied to other current matrix elements where the underlying physics is more 
involved. We plan to study those questions in future publications. 
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